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ABSTRACT 

The Chandra image of the central region of the Coma cluster reveals that both its dominant galaxies, NGC 4874 
and NGC 4889, retain the central parts of their X-ray gas coronae. The interstellar gas with a temperature of 1- 
2 keV is confined by the hot intergalactic medium of the Coma cluster into compact clouds (only 3 kpc in radius) 
containing 10** Mq of gas. The physical state of the gas in these clouds appears to be determined by a delicate 
balance between radiative cooling and suppressed (by a factor of 30-100) heat conduction through the interface 
between these clouds and the hot cluster gas. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Coma) — magnetic fields — 
X-rays: galaxies — galaxies: individual (NGC 4874) — galaxies: individual (NGC 4889) 



\. INTRODUCTION 

ROSAT observations showed that the Coma cluster, once 

■ thought to be an archetype of a relaxed system, contains sub- 
structures on a wide range of linear scales (Briel, Henry & 
Bohringer 1992, White, Briel & Henry 1993, Vikhlinin, For- 
man & Jones 1997). The detected structures include 100 kpc- 
scale X-ray enhancements around the two dominant cluster 
galaxies, NGC 4874 and NGC 4889 (Vikhlinin, Forman, & 
Jones 1994). Recent XMM observations confirmed these find- 

, ings and also revealed unresolved X-ray sources coincident 
' with these galaxies (Arnaud et al. 2001, Briel et al. 2001). 
' In this Letter, we resolve these sources with Chandra. Both 

sources are extended arcsec-scale (« 3 kpc) remnants of galac- 
' tic X-ray coronae with temperatures 1-2 keV, compressed by 
, the hot, 9 keV, gas of the Coma cluster. Their very existence 

poses interesting physical questions. 

The problem of thermal evaporation of cold clouds embed- 
, ded in a hot medium has been considered in a number of papers 

■ (Cowie & McKee 1977, McKee & Cowie 1977, Balbus & Mc- 
; Kee 1982, McKee & Begelman 1990, Bandiera & Chen 1994 
1 among others). This problem is relevant for the survival of cold 

clouds in multi-phase cluster cooling flows (e.g., Bohringer & 
Fabian 1989, Fabian, Canizares, & Bohringer 1994). The gen- 
eral consensus emerging from these studies, as summarized by 
Fabian (1994), is that heat conduction must be suppressed by 
one to two orders of magnitude below the Spitzer value, oth- 
erwise the radiative cooling of the gas will be overcome by 
heat conduction. However, the direct evidence for suppres- 
sion of heat conductivity in the intracluster medium (ICM) is 
scarce (one of a few examples can be found in Bechtold et al. 
1983). The Chandra observations of the interstellar medium in 
the Coma galaxies presented here allow a direct study of cold 
clouds embedded in a hot intracluster medium. 

We use Hq = 50 km s^^ kpc^', which implies a linear scale 
of 0.68 kpc per arcsec at the Coma distance {d = 140 Mpc). 

2. CHANDRA OBSERVATIONS 

The Coma cluster was observed with ACIS in the fall of 1999 
in a series of 6 pointings with individual exposures of 10 ksec. 
The pointing position was kept the same during these observa- 
tions, with NGC 4874 and NGC 4874 located at off-axis dis- 



tances of 2.8' and 4.0', respectively. The detector was moved 
in the focal plane between the pointings, so that the galaxies 
were observed by both ACIS-I (4 times) and ACIS-S (2 times). 
There were no background flares in the Coma observations. In 
one pointing, NGC 4874 was located in the gap between the 
ACIS-I chips. The total exposure was 48 ksec for NGC 4874 
and 56 ksec for NGC 4889. 

The central part of the composite ACIS-I image is shown 
in Fig. 1. Notice the faint arcmin-scale X-ray enhancements 
around NGC 4874 and 4889, previously detected by ROSAT. 
At the center of these extended structures there are bright 
(500 counts) compact sources coincident with the optical po- 
sitions of the galaxies. The observed 0.5-2 keV flux is 
3.9±0.1 X lO^'^^ergs-^ cm^^ for NGC 4874 and3.2±0.1 x 
lO^''* erg s^' cm^^ for NGC 4889, corresponding to luminosi- 
ties in the same band of 9.1 x 10^" and 7.6 x 10'*'' erg s^*. 

The enlarged X-ray and optical images of NGC 4874 and 
4889 are shown in Fig. 2. The X-ray sources associated with 
the galaxies are extended, but unlike most ellipticals, the X-ray 
emission is much more compact than the distribution of optical 
light. Below we discuss these interesting objects in detail. 

3. SPECTRAL ANALYSIS 

The spectra of the NGC 4874 and 4889 sources were ex- 
tracted within a radius of 7" (which encompasses most of the 
observed X-ray flux). Since we are interested in the intrinsic 
properties of the galaxies, we collected the background spectra 
in the 15"-30" annuli around them, so it includes the true back- 
ground as well as the emission of the Coma cluster The source, 
background spectra, and the response files were generated in- 
dividually for each pointing. We then co-added the observed 
spectra and averaged the response files for each source. 

In the spectra for both galaxies (Fig. 2), one immediately no- 
tices strong emission lines near 1 keV — a key signature of a 
1 keV-temperature plasma. Indeed, the spectra are well fit by 
the Mewe-Kaastra-Liedahl model in XSPEC with the absorp- 
tion fixed at the Galactic value A^^ = 9 x lO'^ cm^^. There- 
fore, the observed X-ray emission is most likely produced by 
warm intersteUar gas. The X-ray emission cannot be due to 
the integrated emission of stars because its spatial distribution 
does not follow the optical light (Fig. 2 and 4). Any significant 
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contribution by low-mass X-ray binaries (LMXB) can be ruled 
out because it also should follow the starlight and have a very 
different (hard) spectrum (Sarazin, Irwin, & Bregman 2001). 
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Fig. 1. — The central region of the composite exposure-corrected ACIS-I 
image of Coma in the 0.5-2 keV energy band. The narrow band running from 
bottom-left to top-right is the chip gap; it disappears in the smoothed image. 

The best-fit temperatures and metal abundances of these 
sources are listed in Table 1 (the uncertainties hereafter are at 
the 68% confidence level for one interesting parameter). The 
average temperature is 1.0 keV for NGC4874 and 1.8 keV 
for NGC4889. The average abundance of heavy elements is 
close to the Solar value in both galaxies. This is significantly 
higher than the metal abundance in the surrounding Coma gas 
(0.25 ±0.02, Arnaud et al. 2001). 

Chandra's angular resolution permits spectral fitting in sev- 
eral concentric annuli in NGC4874. This analysis reveals a 
marginally significant (at the 99% confidence level) decrease 
of the gas temperature towards the center of the galaxy (Fig. 3 
and Table 1). NGC4889 does not show any non-isothermality, 
but the statistical errors are much larger. 

Table 1 
Spectral results 

Object fcT(keV) a (Solar) x^/^-o-f 

NGC4874 1.00 ± 0.04 0.79+g|^ 23.5/16 

NGC4874 (0"-1.5") 0.85±0.06 0.79* 16.5/17 

NGC4874(1.5'M") 1.11 ±0.08 0.79* 14.1/17 

NGC4889 l-82lo:io i-^l^ 17.9/16 

* — fixed. 

The errorbars are at the 68% confidence level for one interesting parameter. 
The 90% lower limit on the average metal abundance is 0.43 and 0.81 for 

NGC4874 and NGC4889, respectively. 

We also checked for a possible fall in the projected tempera- 
ture of the main cluster emission on an arcmin scale around the 
galaxies, as observed by XMM around NGC 4874 (Arnaud et 
al. 2001). We find no indication for cool gas beyond the small 
radius of the galactic sources. In the 10"-1' annuli centered on 
the galaxies, we measure T = 8.9 ± 0.4 and 9.2 ± 0.4 keV for 
NGC 4874 and 4889, respectively; these temperatures are con- 



sistent with Chandra and XMM temperatures farther away, in 
the l'-3' annuli around the galaxies. 

4. SPATIAL DISTRIBUTION OF THE INTERSTELLAR GAS 

The X-ray surface brightness profile of NGC 4874 is shown 
in Fig. 4. A comparison with the point source profile (dot- 
ted line) clearly shows that the observed X-ray emission is ex- 
tended. We fit the observed profiles with the /3-model, ng = 
nQ(l + r^ /rl)~^^/^, truncated at some radius, r^m, which is also 
a free parameter. The truncation of the gas distribution at 
is required because the ISM is embedded in the high-pressure 
cluster gas and, obviously, cannot extend beyond the radius at 
which the internal and external pressures are equal. The distri- 
bution of volume emissivity was computed assuming a constant 
temperature for NGC 4889 and the observed temperature gradi- 
ent for NGC 4874 (Fig. 3). 

We had to account for the Chandra PSF in fitting the sur- 
face brightness profiles because the source size is comparable 
to that of the PSF. The PSF model was obtained from the library 
of raytracing simulations; the PSF model is consistent with the 
image of the point source located « 1' East of NGC 4874. 



Table 2 

Gas DENSITY DISTRIBUTION PARAMETERS 



Object 




r^, kpc 


i3 


'■cut' kpc 


NGC 4874 


0.18 ±0.02 


0.8 ±0.3 


0.5* 


2.7 ±0.4 


NGC 4874, no cut .. 


0.14 ±0.02 


2.4 ±0.7 






NGC 4889 


0.09 ±0.01 


2.0 ±0.7 


0.5* 


3.1±0.3 


NGC4889,nocut .. 


0.10 ±0.01 


3.3±0.7 


1 7+°° 

^•'-0.5 





* — fixed. 



The results of the surface brightness profile modeling are 
summarized in Table 2. We were unable to obtain meaningful 
constraints on j3 because of the parameter correlation and finite 
PSF. Generally, all j3 > 0.4 are allowed. The cutoff and core- 
radii are well-constrained for j3 < 1. For larger j3, the cutoff 
radius becomes much larger than the core-radius and hence has 
almost no effect on the observed profile. For a pure /3-model 
without a cutoff, j3 must be in excess of 1.2 — 1.3, mimicking 
the cutoff at Ir^ — 2rc, which is consistent with the best-fit r^.^. 
Of course, different models correspond to very similar absolute 
gas density distributions in the radial range of interest. 

5. DISCUSSION 

5.1. Pressure Equlibrium of ISM and Cluster Gas 

It is interesting to compare the pressure of the ISM at the 
best-fit cutoff radius with the external pressure of the cluster 
gas. Even though the three-dimensional location of the galaxies 
within the cluster is uncertain, the ICM density in their vicin- 
ity can be determined rather accurately because each galaxy 
has an associated arcmin-scale X-ray enhancement (Vikhhnin 
et al. 1994): n-^^^ = 2.9 x 10"^ cm'^ around NGC 4874 and 
5.7 X 10"^ cm"^ around NGC 4889. The temperature of the 
cluster gas around both galaxies is T^^^ = 9 keV (see above). 
Remarkably, we find that for the observed ISM parameters, the 
pressure ratio at r = r^.^ is p^^^/ p^^^ = 1 . 1 ± 0.3 for NGC 4874 
and 1 .2 ± 0.3 for NGC 4889. Therefore, the pressure of the ISM 
in both galaxies equals the external pressure at a radius where 
the Chandra image independently indicates an ISM density cut- 
off. This is a clear indication that the distribution of the ISM is 
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Fig. 2. — Images and spectra of NGC4889 (a, b, c) and NGC4874 (d, e,f). The X-ray (0.5-2 keV) and optical images have the same angular size, 2' x 2'. Solid 
histrograms is panels c and/ show the best-fit model with parameters from Table 1. 



affected by the external pressure. The derived central ISM den- 
sities are higher by a factor of 2-10 than those normally found 
in isolated elliptical galaxies (n^ ~ 0.01 — 0.06 cm~^; Forman, 
Jones & Tucker 1985), which also suggests compression by an 
external pressure. 

5.2. Origin and Lifetime Scales for the ISM 

Stellar ejecta is the most likely source of interstellar gas since 
the metal abundance in the ISM of both galaxies is signifi- 
cantly higher than that in the Coma ICM. The observed gas 
mass (1.1 ±0.1 X 10** M,. and 1.6±0.3 x 10** Mg inNGC4874 
and NGC 4889, respectively) would be ejected by stars over the 
time = Mgas/ (M^L^pj), where is the stellar mass loss rate 
per unit optical luminosity. Using the generally adapted value 
= 1.5 X 10-"M. yr-> (Faber & Gallagher 1976) and 
the optical luminosily within (Faber et al. 1997), we find 

« 8 X 10** yr in both galaxies. Therefore, stellar ejecta easily 
produce the observed ISM over the galaxy Ufetime. The ISM 
is subject to radiative cooling, thermal evaporation due to the 
surrounding hot gas, heating by stellar winds and supemovae, 
and possibly to ram pressure stripping. 

Ram pressure stripping. The observed X-ray centroids are 
within 0.5 kpc of the optical centers (limited by the Chandra as- 
pect accuracy). This and the azimuthally symmetric X-ray mor- 
phology probably exclude any significant ongoing stripping. 
Stripping will be unimportant if the galaxies move together with 
their surrounding 100 kpc scale gas density enhancements. 



Radiative cooling. The average cooling time of the ISM is 



= 3Mi,askT /IfxinpLx. For the observed ISM parameters, 
- 1 .1 X 10^ yr for NGC 4874 and 3.7 x 10^ yr for 



'cool " -'^"gas 

we find t 



NGC 4889, with a 10 - 20% uncertainty. 

Stellar winds and supemovae. After 10^ years past the burst 
of star formation, the average temperature of ejecta from nor- 
mal stars and supemovae is ~ 0.6 keV (David, Forman, & Jones 
1990). Therefore, stellar mass loss cannot be a source of heat- 
ing for the ISM with T = 1 - 1 .8 keV. 

Thermal evaporation. Following Cowie & McKee (1977), 
one can show that classical heat conduction through the ISM- 
ICM boundaries of our galaxies should be saturated. From 
their equation (64), we find the saturated evaporation timescales 
fevap « 3 X 10*' yr for NGC 4874 and « 2 x 10^ for NGC 4889, 
with a 20% formal statistical uncertainty. 

The timescales for radiative cooling and classical evaporation 
are both very short and the ISM would be destroyed quickly by 
either mechanism, unless the radiative cooling and the heat in- 
flux that drives evaporation balance each other. This would be 
achieved if the conductivity through the ICM-ISM interface is 



a factor of feool/^evap 



30 — 100 below the classical value. A 
plausible reason for such a suppression is a disjoint magnetic 
structure of the ISM and ICM. 

5.3. Thermal Balance in NGC 4874 

If heat influx from the ICM does balance radiative cool- 
ing of the ISM, one expects that at each radius, the X-ray 
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Fig. 3. — Temperature profile of NGC4874. Metal abundance and 
absorption were fixed at their average values. The statistical signifi- 
cance of the temperature difference in the radial ranges 0"-1.5" and 
1.5"^" is 99%. The temperature gradient can be modeled as T/keV = 
0.87 + 0.5 lg(r/arcsec) 

volume emissivity of the ISM equals the divergence of the 
heat flux (i.e., the local heating rate). Interestingly, this can 
be verified for NGC4874, where we observe a temperature 
gradient and therefore can estimate the heat flux. Assum- 
ing that conductivity in the ISM is classical, the heat flux is 
<7 = — 6 X IQ^'' 1'^ dT / dr erg cm^^ s^^ where T is the plasma 
temperature in degrees (McKee & Cowie 1977). For the ob- 
served radial dependence of temperature in NGC 4874 which 
can be modeled as T{r) = 1.0 x 10^ [H-0.71g(r/0.68kpc)] 
(Fig. 3), the divergence of q is indeed within a factor of 1.5 
of the observed X-ray volume emissivity. 

For NGC 4889, the uncertainty in the temperature profile is 
too large to compute the heat flux. However, we can say that its 
ISM temperature is too high to be plausibly maintained by in- 
ternal heat sources such as steUar mass loss and supemovae, so 
there should be external heating, most Ukely due to conduction 
through the ISM boundary. 

6. SUMMARY 

Chandra observations show that the dominant Coma galaxies 
retain a fraction of their X-ray coronae in the form of compact 
(^ 3 kpc in size) clouds with temperatures of 1.0 and 1.8 keV 



Fig. 4. — X-ray surface brightness profile of NGC 4874 in the 0.5-2 keV 
band. Dotted line shows the expected radial profile for a point-like source; 
it is consistent with the observed profile of the source located 1' to the 
East of NGC 4874 (Fig. 2). The solid line shows a j3-model fit. The long- 
dashed line shows a /3 -model fit truncated at r = 4". The fits are almost 
indistinguishable because of the PSF smearing. The dashed line shows 
the HST optical surface brightness profile multiplied by the observed ra- 
tio ^LMXE/^opt (Sarazin et al. 2000) and therefore should represent the 
LMXB contribution to the observed X-ray flux. 



which are in the pressure equilibrium with the hot intracluster 
gas. Time scales for radiative cooling and classical evapora- 
tion are very short but the ISM observed in two galaxies is un- 
likely to be a transient phenomenon. We propose that cooling 
in the ISM is balanced by heat influx through the cluster/ISM 
interface. For this, thermal conductivity through the interface 
should be 30-100 times below the classical value, perhaps be- 
cause of the disjoint magnetic field structure. Furthermore, 
in NGC4874, the observed radial temperature gradient implies 
that if the conductivity inside the ISM cloud is classical, then 
heat from outside the cloud is distributed through the cloud vol- 
ume in such a way that it precisely balances radiative cooling 
at each radius. This leaves httle or no energy for the ISM evap- 
oration, preserving the ISM clouds for a sufficiently long time. 

We thank E. Churazov for useful discussions. The financial 
support for this work was provided by NASA grant NAG5-92 17 
and contract NAS8-39073. 
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